ICHEP 94 Ref. gls0084 
Submitted to Pa 19 
PI 11 



OCHA-PP-43 
NDA-FP-15 
June (1994) 



Swimming of Microorganism 
and the String- and Membrane- like Algebra 



Masako KAWAMURA, Akio SUGAMOTO 

Department of Physics, Faculty of Science 

Ochanomizu University 
1-1 Otsuka 2, Bunkyo-ku, Tokyo 112, Japan 



and 

Shin'ichi NOJIRI 

Department of Mathematics and Physics 
National Defence Academy 
Yokosuka, 239, Japan 



Abstract 

Swimming of microorganisms is further developed from a viewpoint of strings and 
membranes swimming in the incompressible fluid of low Reynolds number. In our 
previous paper the flagellated motion was analyzed in two dimensional fluid, by using 
the method developed in the ciliated motion with the Joukowski transformation. This 
method is further reflned by incorporating the inertia term of fluid as the perturba- 
tion. Understanding of the algebra controlling the deformation of microorganisms 
in the fluid is further developed, obtaining the central extension of the algebra with 
the help of the recent progress on the VFi+oo algebra. Our previous suggestion on 
the usefulness of the A^-point string- and membrane-like amplitudes for studying the 
collective swimming motion of — 1 microorganisms is also examined. 



Introduction 



In our previous paper we have studied swimming of microorganisms viewed from string and mem- 
brane theories. Our target is stiU the understanding of the foUowing ways of microorganisms' swimming 
|^:There exist only three different universahty classes of the swimming ways of microorganisms; (1) 
Swimming with cilia is adopted by the spherical organisms with the length scale of 20 2 x 10'*/im, 
an example of which is Paramecium; (2) the smaller microorganisms with the size of 1 '--^ 50/im swim 
with flagella, an example of which is the sperm; (3) the bacteria with the size of 0.2 ^ 5/im swim with 
bacterial flagella the motion of which resembles the screwing of the wine-opener. Why is it possible that 
such simple classification be realized in the swimming problem of microorganisms? 

Recapitulating our previous workjl], we will comment on a few newly obtained results. Starting from 
the work by Shapere and Wilczek||^, we have developed the ciliated and flagellated motions and studied 
the microorganisms' swimming from the viewpoint of algebraic structure existing in the deformation 
operation of the microorganisms in the low Reynolds liquid, as well as from the swimming dynamics of 
a group of microorganisms, which may be connected to the A^-point correlation function of string and 
membrane theories. 

The microorganisms with the length scale L ^ 1, the Reynolds number R satisfies i? ^ 1, so that 
the hydrodynamics in this case leads to the following equations of motion for the incompressible fluid: 

V-v-0, (1) 

and 

Av = -Vp, (2) 

or equivalently 

A(V X v) = 0, (3) 

where p is the pressure and v(a:) is the velocity field of the fluid. The surface of a microorganism 
swimming in the D = 3 dimensional fluid of the real world forms a closed membrane at a fixed times t, 
the position of which can be parametrized by introducing (Z? — 1) parameters (i = 1, • • • , I? — 1) as 
— X'^{t; f ^, • • • , It is sometimes instructive to consider the D — 2 dimensional fluid. Then, 

the surface of a ciliate (flagellate) becomes a closed (open) string and its position can be described by a 
complex number 

Z = x^ +ix'^ = Z{t;e), 

with — TT < < t:. In the sticky fluid of /i ^ 0, there is no slipping between the surface of a microorganisim 
and the fluid, namely, we have the matching condition 

v(x = X(t;0)=X(i;e), (4) 

or in the general coordinate system with the metric tensor g^'^{x) 

g'"'{x)v^{x%^xit-i)^X%t;0. (5) 



A The ciliated motion in the two dimension fluid 

The ciliated motion can be viewed as a small but time-dependent deformation of a unit circle in a properly 
chosen scale, 

Z{t,e) = s + a{t,s), (6) 

where s = e'^ and a(t, s) is arbitrary temporally periodic function with period T satisfying |a(t, s)| <C 1 
with — TT < < n. The complex representation of the velocity vector can be denoted as 

2vg{z,z) = (vi +iv2){z,z) (7) 
2vz{z,z) = (vi - iv2){z,z). (8) 
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By estimating the translational and rotational flows at spacial infinity caused by the deformation of the 
cilia, we have obtained 0{a^) expression of the net translationally swimming velocity w^''"*^ of the ciliated 
microorganism as follows: 

r, (cilia) . /,\ 

2v}p = —ao{t) 

+ ^ 7i(d„Q!_„+i - d„a„_i - d„a-„+3) - ^ ndnan-i, (9) 

n<l Ti>l 

where an(t) is defined by a(t, s) — '^n=-oo'^"i^)^"- other hand, the net angular momentum 

^(ciha) gg^jjjgjj -\yy |-]^g ml cr o Or g aul sm from the fluid becomes 

- ^ n(d„a_„+2 - anan - d„Q;_„+2) + ^ ndnO^T > • (10) 

n<l n>l I 

The net translation and rotation resulted after the period T come from O(a^) terms since the 0{a) terms 
cancel after the time integration over the period. 



B The flagellated motion in two dimensional fluid 

Microorganisms swimming using a single flagellum can be viewed as an open string with two endpoints, 
H and T, where H and T represent the head and the tail-end of a flagellum, respectively. Our discussion 
will be given by assuming that the distance between H and T is time-independent and is chosen to be 4 
in a proper length scale. This assumption can be shown to be valid for the flagellated motion by small 
deformations in the incompressible fluid. Then, at any time i, we can take a complex plane of z, where 
H and T are fixed on z = 2 and —2, respectively. This coordinate system z can be viewed as that of 
the space of standard shapes of Shapere and Wilczek. Time dependent, but small deformation of the 
flagellate can be parametrized as 

z{t,e) ^2{cose + ismea{t,e)), (ii) 

where the small deformation a{t, 9) can be taken to be a real numbeij^ satisfying 

a{t,e) ^ -a{t,-9). (12) 

Here, we parametrize the position of the flagellum twice, starting from the endpoint T at = — tt, coming 
to the head H at = 0, and returning to T again at 9 — tt. Motion of the two branches corresponding to 
— TT < 9 < and tt > 9 > should move coincidentally, which requires the condition ([l^). The Joukowski 
transformation z = z{w) = w -\- , separates the two coincident branches in the z plane to form 
lower and upper parts of a unit circle in the w plane, outside domain of which we are able to study the 
swimming problem of the flagellate in a quite similar fashion to that of the ciliate. The parametrization 
of our microorganism in the w plane corresponding to Eq.(|ll|) is now 



W{t, 9) = e'\l + ait, 9)) + Oia^). (13) 



Using the mode expansion satisfying Eq.(12) 



a{t, 9) = an (t) sin n9, (14) 



^When a is taken to be a complex number, the length of the flagellum is locally changeable at 0{a). 
For such an elastic flagellum, we have similar results to that of the ciliated motion. In case of real a, its 
length is locally preserved at 0(a), giving a non-elastic flagellum, which is the more realistic one. 
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we are able to determine the net swimming velocity tij.(f^'^sciia) gained by the flagellate motion of microor- 
ganisms: 

2^(fiagcUa) ^ _ ^ mamam+1 + ^ mamdm_i, (15) 

m>l m>2 

On the other hand, the angular momentum v fj^^^^^^^^^^ is given by 

After the time integration over the period T, ■uj^^^'^sciia) vanishes since in our first order approximation, 
the length of the flagellum is fixed in the incompressible fluid. Therefore the second order approximation 
is necessary for the non- vanishing v^(*i^s°i^'*). 



C The selection rules and the symmetry of microor- 
ganisms' swimming 

Even though the resuhs in Eqs.(§), and (H) are obtained perturbatively, we are able to read 

from them the characteristics of the microorganisms' swimming; In order for the ciliates to swim or 
rotate, they need the coexistence of the two different Fourier modes of ni, 712. The selection rules for the 
allowed (?^l,?^2) conbinations are 



i) \ni + ~ 3 or |ni ± ri2| = 1 
a) ni = n2 or \ni + —2 



for the cililated translation 
for the ciliated rotation, 



The corresponding selection rules for the flagellate motion are 



in) 



\ni 



"2! = 1 



for the flagellate translation. 



(17) 
(18) 



(19) 



where the Fourier modes are sin n9 in this case. 

Viewing these selection rules, we are tempted to elucidate the algebraic structure possibly existing in 
the background of the swimming mechanism. It is similar to the Virasoro algebra, but is different from 
it. For such a purpose, introduction of the "action" will be convenient. The "action" S reproducing the 
classical equations of motion of the swimming of N microorganisms in the incompressible liquid with low 
Reynolds number may be given by 

Sn = Y.h^ I d''-'iii)Pl:\t-^f^^) fe)(t;e(.))-«n^w(i;ew)) 



j=l 



where the velocity field uj^^ix) is given by 



L M 4 



(20) 



(21) 



whose D = 2 expression is Wz,z(^, z) = dzVg — dgVz- 

We have introduced the parameter a' so as to make S dimensionless, where a' has the dimension 
of (Length)^ (Time) The reason why we have used the notation a', familiar in the string theories 
to describe the Regge slope, will be understood later. The Lagrange multipher fields of P^i^\t;^(i)) 
[i — 1, - ■ ■ , N, jjL = 1, • • • , D) guarantee the matching condition of (||) or (^ for i-ih. microorganism at 
any time, and the pressure p{x) is also such multiplier giving the incompressibility given in Eq.(|l|). The 
field equation (H) can be easily reproduced. In the action (po|), time t appears only in the first term of 
representing the matching conditions, that is, the time evolution is triggerd only by the self-motion of the 
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microorganisms, of which influence spreads instantaneously over the whole space and causes the change 
of the fluid velocity there. Because of the lacking of the kinetic term, we may call S as the "action" . The 
additional metric contribution such as yg{x) is only relevant for the curved space, an example of which 
has appeared in the flagellate swimming on the w plane. The later discussion is given for the flat metric. 
Now, we will define the following local transformation at a fixed time t: 

= AnX(,)(t;ew)) (22) 

SP^'^ = (23) 
Sv^ix) = A^(a;) (24) 
Sp{x) = k{x), (25) 

where we have assumed that the transformation parameters A''(a;) and k{x) are restricted by the equations 
of motion, 

df^X^ix) ^ 0, and a^A^.(x) = -d,K{x), (26) 



where A^^ = d^X^, — d^Xf^ is the vorticity for A^. Meaning of the transformations (|22|)^(|25D are quite 
simple; the deformation of the microorganisms (|2^) triggers the increase of the velocity field ( p^ and of 
pressure (|5|) so that they can be consistent with the incompressible fluid dynamics of the low Reynolds 
number. It is also important to note that the succession of these time-independent transformations result 
in the time evolution of our problem. Therefore, the transformations (|2^)~(p5|) resemble the ordinary 
canonical transformation generated by the Hamiltonian. 

The transformation also generates the deformation of the shapes of microorganism (|2^ ) . The generator 
of this deformation per unit time can be written as 



Vl-y'^"-^eA''(X(t;e))^^^, (27) 
which gives the volume (area for D — 2) preserving diffeomorphisms owing solely to the incompressibility 



condition in Eq.(|26|). The second condition in Eg. ([2 61) adds the further restriction on L]^\'^y. By the 
help of the stream function ij\[x), the incompressibility condition is automatically satisfied through 
— ei^^^duOx for D = 3, so that the Eq.(|2^) becomes the constraint on the stream function 



A(Ag^, -a^a,)a''(a;) -0, (28) 

or in D = 2, in terms of the only non- vanishing component a ~ 

A2cr(x) = 0. (29) 

This constraint has been already solved generally in I? = 3 and D = 2 fluid. 
The deformation operator of the fluid corresponding to ( pTj ) 

becomes in = 2 as 

= 2{XA + X,d,) 

= 2{d,ad,~d,ad,), (31) 

where the stream function a{z,z) contains zz'^, z'^z, Inz, Inz, and zlnzz terms (fc : integer). Then, L„ 
can be considered as a Liouville operator of a dynamical system moving in the phase space of {z,z), 
having —a{z, z) as its Hamiltonian. Invariance of the phase volume during the temporal evolution of the 
dynamical system shows that Lu is really the area preserving diffeomorphism. The commutation relation 
is simple, namely 

[Lai, ^0-2] = ^-^{cri: o-2}i (32) 
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where the {ci, (72} is the Poisson bracket defined by 



{cti, cr2} = dzaidg<72 - dgaidz(T2. (33) 

Recently we have tried to obtain a closed algebra including In^ and Inz in a{z,z). The algebra so 
obtained consists of L(^i^rn,n) = -^z'2™(in22)"j ^^id M = Linz-in?; They satisfy 

+ {mr + np~'lr- nq) i(;+p_i, m+q-i, n+r-i) 

+ c{m - I) 5l+p^o6m+q,o6n+r,0 , (34) 

[M, M] = , (35) 



and 



-^(/,m,n)] = -(^ + "l) L(/_i,m_i,„) - 2nL(i_i^„_i_„_i) 



+ 2'^'^'.0 ^m,0 <5n,0- (36) 

In the above expression we add the central charge c, corresponding to the possible central extension of the 
algebra in which the Jacobi identities are kept to hold and the generators are understood to be properly 
redefined. 

The reason why the Inz or Inz is permitted in the stream function cr(z,z) is that the existing sin- 
gularities at z = can be hidden inside the body of the microorganism itself. Therefore, if we are 
not interested in the circulation flow of the fluid (topological flow) around the microorganism we can 
ignore the logarithmic contribution in a. In that case, the algebra becomes wi+oo'- With the notation 

Tn,m. = L(n, 771,0), we have| 

\I'n,m, Tk,l\ = -{nl - mk)Tn+k-l, 771+1^1 

+c{ra-n)5.,n+k,o5m+i,Q- (37) 

There is a discrepancy between the definition of the Wi+tx; and the above expression (|3^). Corresponding 
to the classical generator z'^D^ of Wi+^a algebra, we define the quantum version of the generator as 

d 

Then, the operator r„,„ = a"+'""2j4^(z""'"D'") satisfy the algebra | 



oz 



in the limit of a ^ for c = 0. 

Recently the representation theory of this Wi+oa algebra is progressing considerably |Q, ||, ^, 0, ^ |^ so 
that we are wishing to apply it to our problem. Successive application of the infinitesimal deformation on 
the microorganism forms a swimming motion. On the other hand a representation of the Wi+oo algebra is 
obtained by a successive operation of raising operators V^^ (n>l,r>l)ona highest weight state |A > 
which is characterized by a set of eigenvalues for Vq (n = 1, 2, 3, . . .). This conceptual correspondence 
between the swimming motion of microorganisms and the representation theory of the algebra, both of 
which are controlled by the area-preserving diffeomorphisms, may be useful to understand algebraically 
the typical pattern of the microorganisms' swimming. 



D Small but Non- Vanishing Reynolds number 

Next, we will study the case in which the Reynolds number R is small but non- vanishing. This has not 
been studied in our previous work [Q. In fact the length size L of the microorganisms ranges from /im 

^ For the ordinary VFi+oo algebra, m > but in our algebra of microorganisms' swimming, m can be 
negative. 
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to mm. Therefore, for the larger microorganisms oi L = O(lmm), R of 0(1) can not be ehminated from 
the beginning. The incompressible fluid with R. =/= has the following equation of motion: 

-{v + (v- V)v} = --Vp + Av, (38) 

or 

e{V X v + V X ((v V)v)} = A(V X v), (39) 

with e = This equaiton reduces to eq.(^ or (|^) for R — Q. We will use e as an expansion 
parameter of incorporating small perturbation from the non-vanishing R. Expansion of v in terms of e, 
V = v(") + v(i) H , leads to 

A(Vxv(")) = (40) 
A(Vxv(i)) = e{V X + V X ((v(") • V)v°)} (41) 

for = 3, but for D = 2 they can be written as 

'id,d,{d,v^^'> - d.v^J"^) = (42) 



TKy.'iYiCy.fiZ 

n<l,m+n>2 n<l,m>n 



-m+n— 1 



(43) 



The matching condition ^ can be understood as 

v(")(x = X(t,0)=X(i,C) (44) 

and 

v(i)(x = X(t,e))-0 (45) 

As an example we will take up the ciliated motion in the D = 2 fluid. Then v'^^ satisfying (^) and 
( ^ ) reads as before 

4°^ (z, z) = v^-^ (z) + w W (z) + iz-^~ z)v'^-\z) (46) 

where 

+ ^ na;;rd„z-"+"+^ + , (47) 

n>l,m>n n<l,m+n>2 I 



n>l,m<n n<Cl.in+n<2 J 

Substitution of v'"' into eg. (^l|) , we can find v*^^) satisfying eq.(^5|) as well as the boundary condition at 
spacial infinity where v'^-* is at least finite. This is carried out by adding properly the arbitrary solution 
v^^^^ satisfying the homogeneous equation A(V x v*^^^ ) = 0. We have found, however, unwanted terms 
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behaving z^z ^, z^z ^, In 2, and Inz in the solution v*-^-*. To ohminate these unwanted terms, we impose 
the following restriction on Q!(s, i); 



dfc + ^ f^{(a_„+fe+ia„ + d_„+fe+id„) 

n<0 



-("n+fc-ian + ari+fe-ia„) - (Q!_„_fe+3d„ + a-„-fe+3Q;„)} 
+ 51 + ttri-fe+ldn) 



n>l 



|^(n-fc-l)(n-fc)d„d_„+fc+i =0, = -2) (49) 



k 



^ n(n + l)d„d_„ = 0, (50) 



n<l 



(/c + l)(d2dfc+i +d_idfe) = (A:<0) (51) 

didfe = (fc < 0) (52) 

~ (cilia) 



Now, we have the final solution of v^^^ +v(^^ This gives the net swimming velocity w^' and the angular 

-(ci 



momentum gained Wp^''"'^ when the microorganism is swimming in the fluid with small but non- vanishing 



Reynolds number; 



-(cilia) . 



k{aka-k+i - dfcafe_i - aka-k+3) - ^ kdkak-i 

k<l k>l 

e i f^ir^ 2(fc-3) ^ 2 . . 

4 I 1^1.^ (fc-l)(fc-2)"'="-'=+^ - 1^"'="-^+^ 

^ 2(3fc + l) . -. A 

(a-fc+id/c - a/c-idfe + afc+idfc - a-fc+ad/c) 

fc<-i 

-(a_2d-i + ald-i)} , (53) 



-(ciha) _ —lT[n{ai{t) 



- ^ fc(dfca_fc+2 ^ dfcQffe - d/ca-fc+2) + ^ kdkOLk 
k<\ k>\ 

^4^™ i X/ ^ _ (d^A;+2dfc - d_fe+2dfc) - (a-id-i + oad-i) 
U->i 

^— j-Y(a-fc+2dfc - Qffcdfc + Ofcdfc - aZFr2dfe) > (54) 



fc<-i 



It is interesting to note that the selection rules (|T^) and (|l^) are also valid in this case of small but 
non-vanishing Reynolds number. 
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E Collective motion of microorganisms and A^-point 
correlation function 



Our previous study in this problem is summarized in the foUowing. Consider the situation where the 
vortices are created and annihilated, so that the probability of having the vortex distribution (field) 
uj^'^{x) is given by 



P[uj'"'{x)] -exp 



{-2^ / d''xVW)l^,A^)^^''i^)], (55) 



where ia' is the external parameter controlling the fluctuation of the vortex distribution. [The a' ^ 
limit corresponds to the classical limit without the fluctuation.] In this situation we should sum over 
all the possible configurations of the velocity fields with Eq.(j5^) as their probability. The probability 
of having N microorganisms with their surfaces located at X(2), • ' ' : ^(iv)i and with their time 

derivatives X(2): ' ' ' tX{N)^ is given by the following A'^-point correlation function (amplitude); 

G'Ar(X(i), • • • ; X(jv), -'^(Ar)) 

xGAr(X(i),p(i); ••• ;X(jv),pW), (56) 

where 

Gjv(X(i),p(^^ ••• ;^W,P''^^) 



cxp 



X pW(t(,);ew)Gf (X(,)(e(,)) -X(,)(e(,))) p^p{tur,^u))] > (57) 

with the Green's function of the transverse waves. Then, we have obtained A^-point correlation function 
for the collective swimming of N microorganisms. It is quite similar to the iV-point function of strings 
for D — 2 case and membranes for D = 3 case. It is also related to the string field theory, since the 
incoming and outcoming strings are not point-like, but the Reggeons. 

Here we will discuss the use of A^-point correlation function (|5^) in the collective swimming motion 
of microorganisms. If Gjv represents the probability of having N microorganisms whose surfaces are 
located at with velocity {i — 1, • • • ,N), then it can be viewed as the probability distribution 
of Xn of the imaginary microorganism N located at spacial infinity under the given data of Xi and 
Xi{i = 1, ■ ■ ■ , N — 1). Following the usual strategy, the counterflow —X^ can be indentified to the 
collective swimming motion of A'^ — 1 microorganisms. Therefore, the averaged collective swimming 
motion over the fluctuation distribution is given by 

— (Xn) — — ^ XnGn{Xi,Xi; • • • ; Xn, Xn)- (58) 

It is also an interesting problem to include the kinetic terms of the microorganisms themselves, or to 
consider the stochastic behavior in the swimming motion of the microorganisms. 



Summary 

In this paper, recaptulating our previous work on the swimming of microorganisms, we have given a 
few newly obtained results: we have obtained the closed algebra and its central extension, which is a 
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generalization of the Wi+oo algebra, controlling the swimming motion of microorganism. For the ciliated 

motion, we have extended our previous result to the larger microorganisms for which the Reynolds number 
is small but non-vanishing. We have developed the perturbation theory with respect to the Reynolds 
number. 
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